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Rotating properties of Bose-Einstein condensates in synthetic magnetic field are studied by numerically solv¬ 
ing the Gross-Pitaevskii equation and compared with condensates confined in the rotating trap. It seems that it is 
more difficult to add large angular momentum to condensates spined up by the synthetic magnetic field than by 
the rotating trap. However, strengthening the repulsive interaction between atoms is an effective and realizable 
route to overcome this problem and can at least generate vortex-lattice-like structures. In addition, the validity 
of the Feynman rule for condensates in synthetic magnetic field is verified. 
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I. INTRODUCTION 

Rotating quantum gas has attracted enormous research in¬ 
terest since it exhibits a number of counter-intuitive phenom¬ 
ena in contrast to the classical gas. One of the striking features 
is that vortices inside are quantized. Such features have been 
extensively explored in the context of superfluid ^He [1] and 
^He [2, 3], and superconductors [4]. The recent achievement 
of atomic Bose-Einstein condensates (BECs) opens up a new 
avenue towards understanding the physics of rotating quantum 
gases [5]. The atomic BEC has advantages in comparison to 
superfluid ^He and ^He, since it is flexibly manipulated and 
can be easily rotated by means of several techniques. 

The first reported vortex in BECs was created by the phase 
engineering technique with laser beams, in which a binary 
mixture of condensates was manipulated and quantized rota¬ 
tion of one component was realized [6]. Soon after, the rotat¬ 
ing frame method succeeded in experiments, where the con¬ 
densate was spined up by the rotating deformation of the con¬ 
fining trap [7]. This approach is similar to the rotating-bucket 
method in experimental studies of superfluid "^He [1]. One 
quantized vortex appeared as the rotating frequency rose to a 
certain critical rotating frequency of the deformed trap [7, 8]. 
With the increasing of rotating frequency, the number of vor¬ 
tices increases larger and larger. Vortex lattices [9-11] are also 
obtained for higher rotating frequency. Eurther increasing ro¬ 
tating frequency will drive the condensates transforms to fast 
rotation regime [12, 13]. 

The angular momentum of a rotating classical fluid is pro¬ 
portional to its rotating frequency. Eor a quantum gas, the an¬ 
gular momentum is characterized by the total number of vor¬ 
tices it carries. The relation between the number of vortices 
and the rotating frequency is described by the well-known 
Eeynman rule [14], which can be expressed as 2'KhNv/m = 
2flA, where fl is rotating frequency and Ny is the vortex 
number within area A. The Eeynman rule was deduced origi¬ 
nally for the superfluid helium in the rotating-bucket. In view 
of the analogy between the rotating-bucket and the rotating 
trap, the Eeynman rule is naturally applicable to the atomic 
BECs in the rotating trap. The validity of the Eeynman rule 
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has been intensively studied both theoretically [15-17] and 
experimentally [10, 18] in recent years. 

The rotating frame approach is subject to some limitations, 
e.g., it is difficult to add optical lattices and the rotation is 
limited by heating, metastability, etc. The synthetic magnetic 
held approach stands out and is expected to break above lim¬ 
itations [19, 20]. This novel approach creates vector potential 
for atoms by dressing them in a space-dependent manner with 
optical held and thus makes neutral atoms behave like charged 
particles in the magnetic held [19]. Vortices were observed 
experimentally in the condensate for synthetic magnetic held 
greater than the critical value [20]. Hydrodynamical behaviors 
of the condensate in synthetic magnetic field and the dynam¬ 
ical instability of vortex nucleation were studied in Ref. [21]. 
Enormous effort was devoted to studying cold atoms subject 
to synthetic magnetic field in the presence of optical lattices. 
Experimentally, the synthetic magnetic held has been engi¬ 
neered in periodic lattices [22]. Theoretically, Berezinskii- 
Kosterlitz-Thouless transition in the two-dimensional lattice 
has been investigated [23, 24] and the Hofstadter butterfly 
physics in the strong field has been discussed [25, 26]. 

Thermodynamic properties of an ideal Bose gas in the syn¬ 
thetic magnetic held have also been studied [27, 28]. It is 
natural to expect that rotating the gas leads to decrease of the 
Bose-Einstein condensation temperature. It shows that Bose- 
Einstein condensation can be more easily suppressed in ro¬ 
tating frame than in synthetic magnetic held, which implies 
that the rotating frame can spin up the atomic gas more effi¬ 
ciently than the synthetic magnetic held. If it is true, it is more 
difficult for the synthetic magnetic held to add large angular 
momentum to the condensates than the rotating frame. There¬ 
fore, it is worthwhile to check this issue carefully and to make 
an elaborate comparison of the two approaches. 

In this paper, we focus on rotating properties of atomic 
BECs in the synthetic magnetic held. The vortex formation is 
investigated by numerically solving the Gross-Pitaevskii (GP) 
equation, with the emphasis on the difference between the re¬ 
sults of the synthetic magnetic held and the rotating frame 
approach. This paper is organized as follows. Section II in¬ 
troduces the Gross-Pitaevskii equations of the BECs in both 
the magnetic held and the rotating frame. Section III presents 
the numerical results and discusses the difference between the 
two approaches. A brief conclusion is given in Section IV. 
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II. THE MODEL 


For a neutral atom of mass m in synthetic magnetic field, 
the effective Hamiltonian is given by 

H = ^{P - A )2 + ^ y2^ ^ ( 1 ) 

where P is the canonical momentum operator and A is the 
synthetic gauge potential. Here the atomic Bose-Einstein con¬ 
densate is confined in the anisotropic harmonic trap, and ojo 
and ujz are the trap frequencies in the x-y plane and the z axis. 
We assume that the synthetic field B = 2niPlez is parallel to 
the z axis and choose the symmetric gauge A = (B x r)/2 = 
m {ft X r), where fl is the Larmor frequency. Therefore, the 
GP equation for a Bose-Einstein condensate can be written as 


ih 


dt 


— {P-Af + V'{v) 


■gW 




( 2 ) 



FIG. 1. Contour plots of the density distributions |'!/’ 2 dP showing a 
steady vortex state in a rotating BECs with the dimensionless contact 
interaction <720 = 100. Upper panels correspond to rotating frame. 
The middle (lower) panels correspond to synthetic magnetic field. 
The number in the panels shows the rotating frequency in units of cuq. 
The last figure in the upper panels corresponds to rotating frequency 
0.99a;o. The field of view in the panels is 8ao x 8ao- 


Eor comparison, we also consider properties of BECs in the 
rotating frame. In this case, the single-particle Hamiltonian is 
given by 


where U'(r) = + y"^) + '^(r,f) is the 

wave function of the condensates, and g = Airh^aslrn is the 
contact interaction strength between atoms with s-wave scat¬ 
tering length Qs- In the present work, we consider pancake¬ 
shaped BECs, that is, the harmonic trapping frequencies sat¬ 
isfy UJz > Wq. 

Eor convenience, Eq. (2) can be expressed with dimension¬ 
less quantities, where the spatial coordinates x, y and z are 
normalized by the characteristic harmonic oscillator length 
oq = with 02 = min {uiqjUJz} = ujq and the time 

t is in units of Then the Eq. (2) is reduced to a dimen¬ 
sionless form 


* di 



A)2 + U(f) + /3|Vi|2 




(3) 




(5) 


Here A ~ miVl x r) and H denotes the rotating frequency 
of the trap which is comparable to the Larmor frequency in 
Eq. (1). Therefore, both of them are referred as the rotating 
frequency in this paper. With the same procedure, we have the 
corresponding 2D GP equation in the dimensionless form 


. d'^2DiS:,y,i) 
‘ di 


+32d|'*/'2dP ‘ip2D{x,y,i). 


f) 


( 6 ) 


where V{y) = ^70 + y'^) + 7 o = wo/w = 1, 

72 = uJzjojQ, P = 47 rAas/ao. Here, dimensionless variables 
are denoted with a tilde. Since uJz ^ wq, the wave function 
can be supposed to be in a variables separation form 

, i) = ^ 2 D(i, y, 

where 


In our calculation, we use the Eourier spectral method to 
solve the nonlinear differential Eqs. (4) and ( 6 ) via the imag¬ 
inary time propagation approach [29]. Although these two 
equations have no fundamental difference, the comparison in¬ 
dicates apparent different physics behaviour. The calculation 
starts with some initial states, then propagates until numerical 
convergence is achieved. 


■010(5) = (72/7r)^/'‘e 


is the ground state of the harmonic oscillator along the z di¬ 
rection. After integrating out the coordinate z, we obtain the 
two-dimensional (2D) GP equation 


.9020(5, 

I -=- 

dt 


i(p_i)2 + i(i2+y2) 


+ff2o|02Op 02o(i,y,0, 


(4) 


where 520 = P represents the effective 2D interac¬ 

tion strength. We note that 7 ^ is just the trap aspect ratio of 
the anisotropic trap. 


III. RESULTS AND DISCUSSIONS 


The Gross-Pitaevskii equation provides a remarkably re¬ 
liable description of the physics of the atomic condensate. 
Properties of atomic BECs in the rotating frame have been in¬ 
tensively investigated based on the GP equation [15-17]. The 
formation of vortices with relation to the rotating frequency, 
the contact interaction between atoms and the trap aspect ra¬ 
tio of traps have been discussed in detail. Hereinafter, corre¬ 
sponding results for the BECs in synthetic magnetic field are 
examined. The emphasis is laid on the difference between the 
two cases. 
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A. Vortex formation with relation to rotating frequency 

The formation of quantum vortices and the dependence of 
the vortex number on the rotating frequency are among the 
most important issues for rotating quantum gases. Different 
from the classical case, no angular momentum can be added to 
the quantum gas until the rotating frequency exceeds a critical 
value, when one quantize vortex begins to be created. Then 
more and more vortices appear as further increases the rotat¬ 
ing frequency. 

Figure 1 shows the two-dimensional atom density at differ¬ 
ent rotating frequencies for BECs with the dimensionless con¬ 
tact interaction p 2 D = 100 in the rotating frame and synthetic 
magnetic field. The critical rotating frequency of creating the 
first vortex inside the BECs in the rotating frame and synthetic 
magnetic field is O.SOOwq 0.388aJo> respectively. The two 
values are very similar and the latter is slightly larger than the 
former. 

However, the two rotating approaches take on remarkable 
difference in high rotating frequency cases. Eor BECs in the 
rotating frame, the vortex number grows quickly with H, es¬ 
pecially when n approaches the trapping frequency ujq. Eor 
example, there are about 10 vortices at H/wq = 0.95 and the 
vortex number amounts to 56 at H = 0.99wo which is just 
the upper limit of the rotating frequency. Meantime, the BEC 
expands significantly. It can be seen that the condensate oc¬ 
cupies almost the whole square region under consideration in 
our calculation when H rises to 0.99wo. The vortex structure 
looks like the Abrikosov vortex lattice. We note that the ob¬ 
tained results may not be sufficiently accurate when H is close 
to ujQ, because the boundary condition we choose in the cal¬ 
culation may not be well met. 

Eor the BECs rotated by the synthetic magnetic field, it is 
obvious that the vortex number grows much more slowly. Af¬ 
ter one vortex is created, the vortex number keeps unchanged 
even if il reaches O.Swq. It seems rather difficult to add more 
vortices. Only 4 vortices emerge at O = 0.99a;o. The Larmor 
frequency can be larger than ujq. Actually, there is no intrinsic 
upper limit value for O and it only depends on the strength 
of the synthetic magnetic field. Eigure 1 shows the results for 
O > Wo- The vortex number rises, but still in a slow way. 
There are about 12 vortices at O = l.bwQ. The number is still 
small. In experiments, about 10 vortices were observed in the 
BECs under the synthetic field [19, 20]. 

Another difference from the rotating frame case is that the 
size of the condensates maintains almost unchanged as O in¬ 
creases. This point is clearly shown in Eig. 1. 

Above results indicate that it is hard to create vortex in 
BECs rotated by the synthetic magnetic field, satisfying the 
indication of thermodynamic calculations [27, 28]. To illus¬ 
trate this point further, we calculate the critical rotating fre¬ 
quencies, Hc/wo, of creating certain number of vortices, Ny, 
for the two rotating cases. The results are shown in Eig. 2. As 
mentioned above, in two cases, the critical rotating frequen¬ 
cies for creating one vortex are very similar. As the vortex 
number increases, the two corresponding critical rotating fre¬ 
quencies become more and more different. Eor example, in 
order to create twelve vortices, the critical rotating frequency 



FIG. 2. Plot of the equilibrium vortex number Ny and rotating fre¬ 
quency of creating certain number vortices He/wo for conventional 
BECs with the dimensionless contact interaction g2D = 100. The 
square points and circle points are numerical results. The black and 
red lines represent synthetic magnetic field and rotating frame, re¬ 
spectively. 




FIG. 3. The number of vortices Ny (black square points) and an¬ 
gular momentum per atom i^/h (red circle points) versus rotating 
frequency H/wq for conventional BECs with the dimensionless con¬ 
tact interaction giD = 100. The upper and bottom panels correspond 
to rotating frame and synthetic magnetic field, respectively. 


in the synthetic field case is about twice as large as that in 
rotating frame case. 

It should be noted that the vortices are admitted into the 
condensate via the dynamically instability [21, 30] and if 
the dynamical process of seeding vortices inside the con¬ 
densate is taken into account the critical rotating frequency 
should be larger than the value obtained within the present 
Gross-Pitaevskii equation approach. Nevertheless, as ana¬ 
lyzed above our treatment can provide a qualitatively reason¬ 
able description of the vortex formation process. 

The Feynman rule for atomic BECs in the rotating trap 
has been studied [10, 15-18]. To proceed, we assess the va¬ 
lidity of the Feynman rule for the BECs in synthetic mag¬ 
netic field. Figure 3 shows the dependence of the vortex 
number Ny and the angular momentum per atom on 
the rotating frequency H/wq. The number Ny is counted 
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from the density distribution and the angular momentum 
per atom is numerically calculated according to the equa¬ 
tion 4 = // ilj*Lz'4’dxdy/ JJ lipl'^dxdy, where = 
—ih{xdy — ydx) is the z-component of the angular momen¬ 
tum operator. Numerical results show that the angular mo¬ 
mentum per atom .fz/fi is about a half of the number of vor¬ 
tices, the Feynman rule is basically met for the two frames. 
The small disagreement between Ny and may be at¬ 
tributed to the inhomogeneous density [15]. 


B. Vortex formation with relation to contact interaction 

According to the investigation on BECs in the rotating 
frame, the interaction between atoms can affect the formation 
of vortices considerably [31]. Vortices can be created more 
easily in the BECs with stronger repulsive interactions. It is 
of interest to study the role of the contact interaction for BECs 
in the synthetic field. Eigure 4 plots the critical rotating fre¬ 
quency flci /wq for the single-vortex state as a function of the 
dimensionless contact interaction gro- The lower critical ro¬ 
tating frequency decreases monotonically with the contact in¬ 
teraction for both rotating approaches, which indicates that the 
interaction results in similar effect on the formation of vortices 
in BECs rotated by the synthetic magnetic field. 
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FIG. 4. The lower critical rotating frequency flci/wo for the single¬ 
vortex state as a function of the dimensionless contact interaction 
g 2 D. The square points and circle points are numerical results. The 
black and red lines represent synthetic magnetic field and rotating 
frame, respectively. 


Eigure 5 illustrates vortex structures at different values of 
92 d- The rotating frequency is relatively large, fl = O.Twq- 
The upper and bottom panels show the atom density for the 
rotating frame and synthetic field cases, respectively. By sim¬ 
ply counting, we find the ratio of vortex numbers between the 
two cases are 5/1, 16/9, 24/17, 39/28, and 46/39 at gxo = 100, 
500, 1000, 2000 and 3000, respectively. Eor both cases, the 
vortex number increases with increasing the interaction, and 
the ratio decreases and tends to 1. These results suggest that 
although it is still relatively difficult to generate vortices in 
the latter case in comparison to the former one, the differ¬ 
ence between two cases becomes small at sufficiently large 
contact interactions. Eor example, vortex numbers are com¬ 
parable at P 2 d = 3000 and lattice-like structures are formed 
in two cases. The vortex lattice hasn’t been observed in ex¬ 
periments by the latter approach [20]. In atomic gases, the 
s-wave scattering length can be tuned easily via the Eeshbach 
resonance [32, 33]. So increasing of contact interactions is an 
effective and realizable route to add more angular momentum 
and generate more vortices in the condensates rotated by the 
synthetic magnetic field. 

As shown in Eig. 5, the condensate size is enlarged appar¬ 
ently with strengthening the interaction. It appears that it is 
easier to add vortices in a condensate of larger size when the 
particle number keeps constant. This insinuate the possible 
reason why it is easy to produce more vortices in the BECs 
spined up by the rotating trap than by the synthetic magnetic 
field. This might be due to the confinement of trap potential. 
According to Eqs. (4) and (6), in the rotating frame, the effec¬ 
tive trapping potential is weakened with the rotating frequency 
n, thus the condensate expands correspondingly. Neverthe¬ 
less, in the synthetic magnetic field case, the trapping poten¬ 
tial is independent on the Larmor frequency. 

At last, we briefly discuss the effect of the trap aspect ra¬ 
tio, 7 z = Wz/woj on the formation of vortices. This issue 
has been discussed for 3-dimensional condensates in the ro¬ 
tating frame and it is indicated that increasing 7 ^ helps nu- 
cleation of the vortex [34]. In the present study the 3D sys¬ 
tem is converted into an effective 2D model and the trap as¬ 
pect ratio 7 z is included in the effective 2D contact interaction 
<?2D = /3\/7z/27r, the effect of increasing 7 ^ is similar to that 
of increasing 520 > as shown in Eig. 5. 


IV. SUMMARY 



FIG. 5. Upper and bottom panels: density profiles for rotating frame 
and synthetic magnetic field at fl = O.Tcvq. From left to right, the 
dimensionless contact interaction is 100, 500, 1000, 2000 and 3000, 
respectively. The field of view in the panels is 8ao x Soq. 


We have investigated rotating properties of the Bose- 
Einstein condensates in synthetic magnetic field. The for¬ 
mation of vortices are calculated by numerically solving the 
Gross-Pitaevskii equation, considering effects of the Larmor 
(rotating) frequency, the interaction between atoms and the 
trap aspect ratio. The obtained results are compared with 
those of the BECs in the rotating frame. We find that the 
vortex number in the condensate spined up by the synthetic 
magnetic field is much smaller than in that by the rotating 
frame given the same rotating frequency, which implies that 
the synthetic magnetic field has less efficiency in rotating the 
condensate. Strengthening the repulsive interaction between 
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atoms or increasing the trap aspect ratio are helpful for cre¬ 
ating more vortices. In particular, when the interaction is 
sufficiently strong, comparable large number of vortices are 
produced by both rotating approaches. Abrikosov lattice-like 
structures can be formed in two cases. Mover, the validity of 
the Feynman rule is checked for the BECs in synthetic mag¬ 
netic field case. 
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